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ABSTRACT 

Anomalous X-ray pulsars (AXPs) and soft gamma repeaters (SGRs) are two small classes of X- 
ray sources strongly suspected to host a magnetar, i.e. an ultra-magnetized neutron star with B k, 
10^'*-10^^ G. Many SGRs/AXPs are known to be variable, and recently the existence of genuinely 
"transient" magnetars was discovered. Here we present a comprehensive study of the pulse profile and 
spectral evolution of the two transient AXPs (TAXPs) XTE J1810-197 and CXOU J164710.2-455216. 
Our analysis was carried out in the framework of the twisted magnetosphere model for magnetar 
emission. Starting from 3D Monte Carlo simulations of the emerging spectrum, we produced a large 
database of synthetic pulse profiles which was fitted to observed lightcurves in different spectral bands 
and at different epochs. This allowed us to derive the physical parameters of the model and their 
evolution with time, together with the geometry of the two sources, i.e. the inclination of the line-of- 
sight and of the magnetic axis with respect to the rotation axis. We then fitted the (phase-averaged) 
spectra of the two TAXPs at different epochs using a model similar to that used to calculate the pulse 
profiles (ntzang in XSPEC) freezing all parameters to the values obtained from the timing analysis, 
and leaving only the normalization free to vary. This provided acceptable fits to XMM-Newton data in 
all the observations we analyzed. Our results support a picture in which a limited portion of the star 
surface close to one of the magnetic poles is heated at the outburst onset. The subsequent evolution 
is driven both by the cooling/varying size of the heated cap and by a progressive untwisting of the 
magnetosphere . 

Subject headings: radiation mechanisms: non-thermal — sources (individual): XTE J1810-197, CXOU 
J164710. 2-455216 — stars: magnetic fields — stars: neutron 



1. INTRODUCTION 

I> ' In recent years an increasing number of high- resolution 
^O ! spectral and timing observations of isolated neutron stars 
^O ■ has become available. Many of these observations con- 
■^ ' cern two peculiar classes of high-energy pulsars, the 
'^ ] Anomalous X-Ray Pulsars (AXPs: 9 objects plus 1 can- 
f^^ , didate) and the Soft Gamma-Ray Repeaters (SGRs: 6 
f^ ' objects) ^. Historically these two classes of sources were 
f^ [ regarded as distinct. While SGRs were first discovered in 
^H , late 1978-early 1979, when SCR 1806-20 and SC R 0526- 
66 exhibited a bright burst of soft 7-rays tiMazets et al.l 
119791 iLaros et a l. 1986), AXPs were observed for the first 
S^ i time in 1981, when .Fahlman. &: Gregorvl (jl981| ) discov- 
Cl ■ ered pulsations in the EINSTEIN source IE 2259-^586. 
Ci ' It was, however, not until the mid '90s that AXPs were 
recognized as a class of "anomalous" pulsars because of 
their lumin osity substantially ex ceeding rotational en- 
ergy losses (jMereghetti. fc Stellaiil995.') . 

Although SGRs were mainly known as emitters of 
short, energetic bursts, they are also persistent X-ray 
sources with properties quite similar to those of AX Ps 
(see e.g. iWoods fc ThomDson|[2006l: lMereghettill2008l for 
reviews). They all are slow X-ray pulsars, with spin pe- 
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riods in a very narrow range (P '^ 
large spin-down rates (P ^ 10^^'^- 

yr, ^ 

compared to those of rotation or accretion powered pul- 
sars {B - 10"-10i5 G > Bqed ^ 4^-4: X lO^^ G). AXPs 
and SGRs have persistent X-ray luminosities Lx ^ 10'^^- 
10'^^ erg s~^. Their spectra in the 0.1-10 keV band are 
relatively soft and can be empirically fitted with a two- 
component model, an absorbed blackbody {kT ^ 0.2-0.6 
keV) plus a power-law (F ~ 2-4). INTEGRAL obser- 
vations revealed the presence of sizeable emission up to 
^ 200 keV, which accounts for up to 50% of the total fiux. 
Hard X-ray spectra are well represented by a power-law, 
which dominates above ~ 20 keV in AXPs. 

The large high-energy output can not be explained 
in terms of rotational energy losses, as in conventional 
models for radio-pulsars, while the lack of stellar com- 
panions argues against accretion. The powering mecha- 
nism of AXPs and SGRs, instead, is believed to reside 
in t he neutron star ultra-strong magnetic field (magne- 
tar; Duncan, fc ThompsonI 119921 [Thompson fc DuncanI 
11993 ). The magnetar scenario appears capable to explain 
the p roperties of both the bursts (Thompson & Duncag 
Il995f ) and the persistent emission (the twisted magne- 
tosphere niodel.JThompson. Lvutikov. fc Kulkarnil [20021 : 
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IZane et all 120091 and references therein; see § [3] for de- 
tails), although no definite model for the hard tails was 
put forward as yet. 

The persistent emission of SGRs and AXPs is now 
known to be variable. AXPs, in particular, display differ- 
ent types of X-ray flux variability: from slow, moderate 



Albano et al. 



flux changes on timescales of months/years, to intense 
outbursts with short rise times (~ 1 day) lasting ^ 1 
year. Some AXPs were found to undergo intense and dra- 
matic SGR-like burst activity on sub-second timescales 
(XTE J1810-197, 4U 0142-1-614, IE 1048.1-5973, CXOU 
J164710.2-455216 and IE 2259-^586). The discovery of 
bursts from AXPs is regarded as further evidence in favor 
of a common nature of AXPs and SGRs. 

The first case of AXP flux variability was observed in 
2002, when IE 2259-1-586 showed an increase in the per- 
sistent flux by a factor ^ 10 with respect to the quiescent 
level, followed by the emission of ~ 8 short bursts wit h 
luminosity Lx ^ lO^^ - 10^8 erg s"! (jKaspi et al.ll2003[ ). 
In early 2003 the 5.54 s AXP XTE J1810-197 was dis- 
covered at a luminosi ty ~ 100 great er than its quiescent 
value (lO^'^ erg s~^: llbrahim et al.l [2004'). Analysis of 
archival data revealed that the outburst started between 
November 2002 and January 2003. 

On September 21st 2006 an outburst was ob- 
served from the AXP CXOU J164710.2-455216 (P = 
10.61 s). The flux level was ^ 300 times higher 
than that measured only 5 days ear lier by XMM- 
Newton (Muno ct al. 20 06b; Campana . fc Israeli 120061 : 
llsrael. fc Campaiial 120061 ) . This, much as in the case 
of XTE J1810-197, indicates that some AXPs are tran- 
sient sources (dubbed Transient AXPs) and may become 
visible only when they enter an active state. Recently 
other AXPs and SGRs showed a series of short bursts 
of soft 7-rays w hich was detected by different satellites 
(jMereghetti et al. 2009 ). 

In this paper we present a comprehensive study of 
the pulse profile and spectral evolution of the TAXPs 
XTE J1810-197 and CXOU J164710.2-455216 through- 
out their outbursts of November 2002 and September 
2006, respectively. By confronting timing data with syn- 
thetic lightcur yes obtained from the twisted magneto- 
sphere model (|Nobili. TuroUa. fc Zand I2008D , we were 
able to estimate how the physical parameters of the 
source (surface temperature and emitting area, electron 
energy, twist angle) evolve in time. The fits of the pulse 
profiles also allowed us to infer the geometry of the two 
systems, i.e. the angles between the magnetic and ro- 
tational axes and the line of sight. Spectral models, ob- 
tained with the parameter values derived from the timing 
analysis, provide acceptable fits to XMM-Newton data. 

2. TRANSIENT AXPS PROPERTIES 

2.1. XTE J1810-107 

The Transient AXP (TAXP) XTE J1810-197 was 
serendipitously discovered in 2003 with the Rossi X-Ray 
Timing Explorer (R XTE) while observing SGR 1806-20 
(jlbrahim et al.l [20041 ). The source was readily identified 
as a X-ray pulsar, and soon after a search in archival 
RXTE data showed that it produced an outburst around 
2002 November, followed by a monotonic decline of the 
X-ray flux. The X-ray pulsar spin period was found 
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to be 5.54 s, with a spin-down rate '^ 10" 
Using the standard expression for magneto-rotational 
losses, the inferred value of the (dipola r) magnetic field 
is B ~ 3 X 10^'' G (jlbrahim et al.ll2004[ ). The source was 
classified as the first transient magnetar. The TAXP 
XTE J1810-197 was then studied with Ch andra and 
XMM-Newton f,Gotthelf et al., ,2004; llsrael et al.. ,2004; 



iGotthelf. fc Halpernll2005ll2007D . in order to monitor its 
evolution in the post-outburst phase. 

By using archival Very Large Array ( VLA) data a tran- 
sient radio emission with a fiux of ~ 4.5 mJy at 1.4 GHz 
was discovere d at the Chandra X-ray position of XTE 
J1810-197 (Hal pern et"aIll2005D . Only later on it was 
discovered that this radio emission was pulsed, highly po- 
larized and with large flux var iability even on very short 
timescales (jCamilo et al.ll2Q06l ). The X-ray and the radio 
pulsations are at the same rotational phase. Since accre- 
tion is expected to quench radio emission, this is further 
evidence against the source being accretion-powered. 

Deep IR observations were performed for this source, 

revealing a weak {Kg = 20.8 mag) counterpart, 

with charac teristics similar to those of other AXPs 

Israel et all (20 04) 



,_. .__ _ _^ . The IR emission is variable 

Rea et al.l 1200411 . but no correlations between the IR 



and X-ray changes were found up to now. The exis- 
tence of_a^_correlation_at IR/radio wavelengths i s uncer- 
tain (i Camilo et aLll2006ll2007l: iTesta et al.ll2008D . 

XTE J1810-197 was observed 9 times by XMM- 
Newton, between September 2003 and September 2007, 
two times every year. The uninterrupted coverage of the 
source during 4 years provides an unique opportunity to 
understand the phenomenology of TAXPs. Earlier obser- 
vations of XTE J1810-197 showed that the source spec- 
trum is well reproduced by a two blackbody model, likely 
indicating that (thermal) emission occurs in two regions 
of the star surface of different size and temperature: a hot 
one [kT = 0.70 keV) and a warm one {kT = 0.30 keV; 
IGotthelf, fc Halpernll2005l ). XMM-Newton observations 
al so showed that the pulse d fraction decreases in time. 

iPerna fc Gotthelj (|2008l ) discussed the post-outburst 
spectral evolution of XTE J1810-197 from 2003 to 2005 
in terms of two blackbody components, one arising from 
a hot spot and the other from a warm concentric ring. 
By varying the area and temperature of the two regions, 
this (geometric) model can reproduce the observed spec- 
tra, account for the decline of the pulsed fraction with 
time and place a strong constrain on the geometry of the 
source, i.e. the angles between the line of sight and the 
hot spot ax is with respect to the sp in axis. 

Recently, iBernardini et al.l ()2009() by re-examining all 
available XMM-Newton data found that inclusion of a 
third spectral component, a blackbody at ^ 0.15 keV, 
improved the fits. When this component is added both 
the area and temperature of the hot component was 
found to monotonically decrease in time, while the warm 
component decreased in area but stayed at constant tem- 
perature. The coolest blackbody, which appeared not 
to change in time, is associated to emission from the 
(large) part of the surface which was not affected by the 
event which triggered the outburst, and is consistent with 
the spectral properties of the source as derived from a 
ROS AT detection before the outburst onset. Finally, an 
interpretation of XTE J1810-197 spectra in terms of a 
resonant con ipton scat t ering model (RCS, see ^ was 
presented bv lReaet aP (|2008D . 

2.2. CXOU J164710.2-455216 

The TAXP CXOU J164710.2-455216 was discovered in 
two Chandra pointings of the young Galactic star cluster 
Wester lund 1 in May /June 2005. The pe riod of the X-ray 
pulsar was found to be P = 10.61 s (|Muno et al.l 12006.) . 
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with a period derivative P = 9 x 10^"'^'^ s s 
I2007f) . The implied magnetic field is B ^ 

In November 2006, an intense burst was detected by 
the Swift Burst Ale r t Telescope (BAT ) in Westerlund 1 
(jKrimni et all 120061: [Muno et al.ll2006D . Its short dura- 
tion (20 ms) suggested that its origin was the candidate 
AXP. However the event was initially attributed to a 
nearby Galactic source, so the AXP was not promptly 
re-observed by Swift. A ToO observation program with 
Swift was started 13 hrs after the burst, displaying a 
persistent flux level 300 times higher than the quiescent 
one. CXOU J164710.2-455216 was observed in radio, 
with the Parkes Telescope. The observation was car- 
ried out a week after the outburst onset, with the intent 
of searching for pulsed emission similar to that of XTE 
J1810. In this case, however, only a (tight) upper limit to 
the radio flux (40 /iJy) was placed (Bur gav et al.l i2006). 

XMM-Newton observations carried out across the out- 
burst onset show a complex pulse profile evolution. Just 
before the event the pulsed f raction was '^ 65 %, while 
soon after it became --11% (jMuno et al.ll2007| ). More- 
over, the pulse profile changed from being single-peaked 
just before the burst, to showing three peaks soon af- 
ter it. CXOU J164710. 2-455216 spectra in the outburst 
state were fitted either with a two blackbody model 
(fcTi - 0.7 keV, kT2 ~ 1.7 keV), or with a black- 
body plus power-law model (kT ^ 0. 65 keV, F ry 2.3 ; 
iMuno et a l. 2007; Isra el et al.l [20071 ) . iRea et all (|2008l ) 
found that a RCS model also provides a good fit to the 
data. 

3. THE MODEL 

It is now widely accepted that AXPs and SGRs 
are magnetars, and that their burst/outburst activ- 
ity, together with the persistent emission, are pow- 
ered by their huge magnetic field. In particular, the 
soft X-ray spectrum (~ 1-10 keV) is believed to orig- 
inate in a "twisted" magnetosphere (V x B ^ 0; 
[Thompson. Lyutikov. fc Kulkarnil 120021 ) . where the cur- 
rents needed to support the field provide a large enough 
optical depth to resonant Compton scattering of ther- 
mal photons emitted by the star surface. Since charges 
are expected to flow along the closed field lines at rela- 
tivistic velocities, photons gain energy in the (resonant) 
scatterings and ultimately fill a hard tail. 

Most studies on sp ectral formation in a 

twisted magnetosphere (iLvutikov. fc Gavriill 120061; 
Fernandez, fc Thompso^ 120071 ; iNobih. TuroUa. fc Zand 



the electron density in the magnetosphere 



2008| ) axe based on the axially symmetric, force- 



free solution for a twisted dipol ar field presented by 
[Thompson. Lyutikov. fc Kulkarni (2002). This corre- 
sponds to a sequence of magnetostatic equilibria which, 
once the polar strength of the magnetic field Bp is fixed, 
depends only on a single parameter: the radial index 
of the magnetic field p {B (x r^^^"^ , < p < 1) or, 
equivalently, the twist angle 
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where Br , Bg and B^ are the spherical components of 
the field, which depend only on r and 6 because of axial 
symmetry. Knowledge of B fixes the current density j = 
(c/47r)V X B, and, if the particle velocity is known, also 
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where e is the electron charge and {(3) is the average 
charge velocity (in units of c). The charge density of the 
space charge-limited flow of ions and electrons moving 
along the closed field lines is orders of magnitude larger 
than the Goldreich- Julian density, noj, associated to the 
charge flow along the open field lines in radio-pulsars. 

In our investig ation we make use of the spec tral mod- 
els presented by [Nobili. TuroUa. fc Zand pOOS . NTZ in 
the fo lowing), who studied radiative transfer in a glob- 
ally twisted magnetosphere by means of a 3D Monte 
Carlo code. Each model is characterized by the mag- 
netospheric twist A(/)jv-S: the electron (constant) bulk 
velocity /3, and the seed photon temperature kT. The 
polar field was fixed at _Bp = 10^^ G. In the applica- 
tions presented by NTZ, it was assumed that the star 
surface emits unpolarized, blackbody radiation and is at 
uniform temperature. Concerning the present investiga- 
tion, the most critical assumption is that of a globally 
twisted magnetosphere, as discussed in some more detail 
in § [5] Taking a constant value for the electrons bulk ve- 
locity is certainly an oversimplification and reflects the 
lack of a detailed model for the magnetospheric currents. 
In a realistic case one would expect that /3 is a function 
of position. However, resonant scattering is possible only 
where the bulk velocity is mildly relativistic. If along a 
flux tube there are large variations of the Lorentz factor, 
only the region where /3 « 0.5 will contribute to scatter- 
ing. Moreover, preliminary calculations of the dynamics 
of charged particles in a twisted force- free magnetosphere 
performed accounting for both electrostatic acceleration 
and Compton drag indicate that /3 is indeed fairly con- 
stant along the central part of a flux tube (Beloborodov, 
private communication) . The assumption of unpolarized 
thermal radiation is not cogent either, since we are not 
interested in the polarization of the escaping radiation 
and the emergent spectrum is quite insensitive to the 
polarization fraction of the seed photons (see, e.g.. Fig. 
4 of NTZ). 

The code works by dividing the stellar surface into 
A^e X Af$ zones of equal area by means of a (cosO, $) 
grid, where 8 is the magnetic colatitude and $ the lon- 
gitude. After a few scatterings photons escape from the 
neutron star magnetosphere and are collected on a spher- 
ical surface (the "sky") which is divided into A^e, x A^*s 
patches, similarly to what is done for the star surface. 
The key point is that the evolution of seed photons from 
each patch is followed separately. This allows us to treat 
an arbitrary surface temperature distribution without 
the need to perform new Monte Carlo runs, by simply 
combining together models from runs with different tem- 
peratures at the post-production level (the geometry is 
shown in Fig. [1]) 

Monte Carlo models are computed (and stored) for the 
simplest geometrical case, in which the spin and the mag- 
netic axes are aligned. As discussed in NTZ, the most 
general situation in which the spin and magnetic axes 
are at an arbitrary angle ^ can be treated at the post- 
production level. If X is the inclination of the line-of- 
sight (LOS) with respect to the star spin axis and a is 
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Fig. 1. — A schematic view of the neutron star, ii and fi are 
the star spin and magnetic axis, respectively. The dashed line 
corresponds to the line-of-sight. The two angles x and ^ are also 
shown. The star surface is divided into three regions: a hot polar 
cap (red), a warm corona (blue) and a colder zone (gray). 

the rotational phase angle, the co-ordinates of the points 
where the LOS intersects the sky can be found in terms 
of ^, X 9'iid a. The pulse profile in any given energy 
band is then obtained by integrating over the selected 
range the energy-dependent counts at these positions as 
the star rotates (see again NTZ for details). In order 
to compare model lightcurves with observations, integra- 
tion over energy is performed by accounting for both in- 
terstellar absorption and the detector response function. 
Actually, the interstellar absorption cross-section a and 
the response function A depend on the photon energy 
at infinity E = E'-y/l — Rs/Rns, where E is the en- 
ergy in the star frame (which is used in the Monte Carlo 
calculation) and Rs is the Schwarzschild radius (we as- 
sume a Schwarzschild space-time and take Rns = 10 km 
and Mns = 1-44 M©). Our model pulse profile in the 
[-El, E2] energy band is then proportional to 



dE exp[-NHcr{E)]A{E)N{a, E) 



(3) 



where Nh is the hydrogen column density and N{a, E) is 
the phase- and energy-depend ent count rate. In the ap- 
plicat ions below we used the 'Morr ison, fc McCammorJ 
()198^ model for interstellar absorption and, since we 
deal with XMM-Newton observations, we adopted the 
EPIC-pn response function. We remark that the Monte 
Carlo spectral calculation is carried out assuming a 
flat space-time (i.e. photons propagate along straight 
lines), so that, apart from the gravitational redshift, 
no allowance is made for general-relativistic effects (see 
iZane. fc Turollal [20061 for a more detailed discussion). 
In particular, no constraints on the star mass and radius 
can be derived in the present case from the c omparison of 
model and observed pulse profiles (see e.g. iLeahv et al.l 



[2008l[2009l) . 

Finally, phase-averaged spectra are computed by sum- 
ming over all phases the energy-dependent counts. Note 
that < ^ < 7r/2, while x is in the range [0,7r] because 
of the asymmetry between the north and south magnetic 
poles introduced by the current flow. 

4. TAXP ANALYSIS 

Our flrst Step in the study of the two TAXPs XTE 
J1810-197 and CXOU J164710.2-455216 was to repro- 
duce the pulse profiles (and their time evolution) within 
the RCS model discussed in §3] The fit to the ob- 
served pulse profiles in different energy bands (total: 
0.5 keV < -E < 10 keV, soft: 0.5 keV < E; < 2 keV, 
hard: 2 keV < _E < 10 keV) provides an estimate of the 
source parameters, including the two geometrical angles 
(, and X- While the twist angle, electron velocity and sur- 
face temperature may vary in the different observations 
(although they must be the same in the different energy 
bands for a given observation), the fits have to produce 
values of ^ and x which are at all epochs compatible with 
one another (to within the errors) in order to be satisfac- 
tory. We then computed the phase-averaged spectra for 
the two sources at the various epochs for the same sets of 
parameters and compared them with the observed ones. 
There are several reasons which led us to choose such an 
approach . The main one is that, as discussed in NTZ 
(see also iZane et al.l [20090 . spectral fitting alone is un- 
able to constrain the two geometrical angles. Moreover, 
lightcurve fitting allows for a better control in the case in 
which the surface thermal map is complex and changes 
in time (see below). 

For the present investigation, a model archive was gen- 
erated beforehand. Each model was computed by evolv- 
ing Npatch = 225, 000 photons for TVe x A^* = 8 x 4 = 32 
surface patches {Ntot — 7,200,000 photons). The pa- 
rameter grids are: 0.1 < kT < 0.9 keV (step 0.05 keV), 
0.1 < ^ < 0.9 (step 0.1) and 0.2 rad < A^^v-s < 1.2 rad 
(step 0.1 rad). Photons are collected on a A^e^ ^ ^*s = 
10 X 10 = 100 angular grid on the sky, and in A^^; — 
50 energy bins, equally spaced in \ogE in the range 
0.1 - 100 keV. 

The analysis proceeds as follows. We first used the 
principal component analysis (PCA) to explore the prop- 
erties of the lightcurves as a population and to select the 
model within the archive that is closest to the observed 
one at a given epoch. This serves as the starting point for 
the pulse profile fitting procedure, which we performed 
by assuming that the whole star surface is at the same 
temperature. The fitting is then repeated first for the 
case in which the surface thermal distribution consists of 
a hot spot and a cooler region, and then by generating a 
new archive with a finer surface gridding, and applying 
it in the case of a surface thermal map consisting of a 
hot spot, a warm corona and a cooler region (see again 
Fig [T|) . Finally, the source parameters derived from the 
lightcurve fitting are used to confront the model and ob- 
served (phase-averaged) spectra. Phase-resolved spectral 
analysis, although feasible in our model and potentially 
important, was not attempted because the decay in flux 
of both sources makes the counting stat istics rather poor 
after the first one/two observations (see lBernardini et al.l 
120091 for more details in the case of XTE J1810-197). 
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Fig. 2. — Principal component representation of the simulated 
lightcurves in our archive (black squares) together with the ob- 
served lightcurves of XTE J1810-197 (red dots) and of CXOU 
J164710. 2-455216 (green dots). All the pulse profiles refer to the 
0.5-10 keV band. 

4.1. PC A 

The principal component analysis is a method of mul- 
tivariate statistics that allows to reduce the number of 
variables Xi needed to describe a data set by introducing 
a new set variables, the principal components (PCs) Zi. 
The PCs are linear combinations of the original variables 
and are such that Zi displays the largest variance, Z2 the 
second largest, and so on. By using the PCs it is possible 
to describe the data set in terms of a limited number of 
variables, which however, carry most of the information 
conta ined in the original sample (see e.g. lZane. fc Turollal 
I2OO6I and references therein). 

Synthetic lightcurves were generated for 32 phases in 
the range [0, 2tt] and for a 9 x9 angular grid, 0° < ^ < 90° 
(step 10°), 0° < X < 180° (step 20°); the archive con- 
tains a total of 136323 models. Once the PCA was ap- 
plied to the lightcurve set, we found that the first three 
PCs (Zi, ^2,^3) accounts for as much as '^ 90% of the 
sample variance. This means that the entire set is sat- 
isfactorily described in terms of just three variables in- 
stead of the original 32 (see IZane. fc Turollal 120061 for 
an interpretation of ^1,^2, ^3). A graphic representa- 
tion of the lightcurves in the archive in terms of the first 
three PCs is shown in Fig. [51 In the same plot we also 
show the PC representation of the pulse profiles of XTE 
J1810-197 and CXOU J164710.2-455216 at the various 
epochs. The points corresponding to observations fall 
within the volume occupied by models and this guaran- 
tees that there is a combination of the parameters for 
which a synthetic pulse profile reproduces the data. The 
PC representation is also used to find the model in the 
archive which is closest to a given observed lightcurve, 
by looking for the minimum of the (squared) Euclidean 

distance Yl,i=ii^i ~ Z°^^Y between the model and the 
observed pulse profile. 

4.2. XTE J1810-197 

We considered eight XMM-Newton observations, cov- 
ering the period September 2003-September 2007 for the 
TAXP prototype XTE J1810-197 (see Table [D for the 
observation log ). Only EPIC-p n data were used, and we 
refer to Bcrnardini et al.l (|2009[ ) , who analyzed the same 



observations, for all details on data extraction and re- 
duction. All the EPIC-pn spectra were rebinned before 
fitting, to have at least 40 counts per bin and prevent 
oversampling the energy resolution by more than a fac- 
tor of three. 



4.2.1. Pulse profiles 

We started our analysis by making the simplest as- 
sumption about the star surface thermal map, a uniform 
distribution at temperature T. Lightcurves were then 
computed in the total, soft and hard energy band for all 
the models in the archive. Once the model closest to 
each observation (and in each band) was found through 
the PCA, we used it as the starting point for a fit per- 
formed using an IDL script based on the minimization 
routine mpcurvef it .pro. Our fitting function has six 
free parameters, because, in addition to the twist angle, 
the temperature, the electron velocity, the angles x and 
1^, we have to include an initial phase to account for the 
indetermination in the position of the pulse peak. Since 
it is not possible to compute "on the fiy" the pulse profile 
for a set of parameters different from those contained in 
the archive, lightcurves during the minimization process 
were obtained from those in the archive using a linear 
interpolation in the parameter space. 

In this way we obtained a fair agreement with the 
observed pulse profiles (x^ < 1.12 in five out of eight 
observations; see Table [j), and the values of the physi- 
cal parameters (A(/)7v-Si /3, T) turn out to be the same 
(to within the errors) for a given epoch among the dif- 
ferent energy bands, as it needs to be. Moreover, the 
evolution of the twist angle and of the surface temper- 
ature follows a trend in which both quantities decrease 
in time as the outburst declines. This is expected if the 
outburst results from a sudden change in the NS mag- 
netic structure, producing both a heating of the star sur- 
face layers and a twisting of the magnetosphere which 
then dies away ("Thomp son. Lvutikov. fc Kulkarnil 120021 : 
[Beloborodov 2009). However, the model is not accept- 
able since we found that the geometrical angles x and 
^ change significantly from one observation to another, 
and even for the same observation in the different en- 
ergy bands (see Fig. [3] where the parameter evolution is 
shown for the three energy bands). The analysis of the 
hard band was not carried out after September 2006, be- 
cause in both the 2007 observations photons with energy 
> 2 keV are only a few and, as a consequence, lightcurves 
are affected by large uncertainties. 

This shortcoming is most probably due to our over- 
simplifying assumption about the NS thermal map. In 
fact, it was shown that in the post-outburst phase the 
surface temperature distributio n of XTE J1810-197 is 
complex and changes in time (jPerna fc Gotthelj 120081 : 
iBernardini et al.l 120091 although a different emission 
model was assumed in these investigations). In order 
to check this, we tried different configurations, start- 
ing with a two-temperature map: a hot cap centered on 
one magnetic pole with the rest of the surface at a con- 
stant, cooler temperature. In this picture both temper- 
atures as well as the emitting areas, are allowed to vary 
in time. While applying the one-temperature model we 
found that both the 2007 observations were reasonably 
well reproduced with a value of the (uniform) surface 
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TABLE 1 
XTE J1810-197 XMM-Newton observations'' 



Label 



OBS ID 



Epoch Exposure time (s) total counts background counts 



Sep03 


0161360301 


2003-09-08 


5199 


60136 


2903 


Sep04 


0164560601 


2004-09-18 


21306 


89082 


1574 


Mar05 


0301270501 


2005-03-18 


24988 


54279 


1760 


Sep05 


0301270401 


2005-09-20 


19787 


21876 


1311 


Mar06 


0301270301 


2006-03-12 


15506 


12296 


1197 


Sep06 


0406800601 


2006-09-24 


38505 


23842 


2974 


Mar07 


0406800701 


2007-03-06 


37296 


21903 


2215 


Sep07 


0504650201 


2007-09-16 


59014 


34386 


4117 



SepOS 



EPIC-pn 

Sep04 MarOS SepOS Mar06 SepOB Mar07 Sep07 



O.B 
? 0.8 
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Fig. 3. — Parameters evolution for XTE J1810-197, uniform sur- 
face temperature; results refer to the total (red dots), soft (blue 
dots) and hard (green dots) energy bands. Parameter errors are 
calculated by the minimization routine mpcurvef it .pro, and are 
at la. Time is computed starting from the September 2003 obser- 
vation. 

temperature ^ 0.15 keV, comparab le to the quiescent 
one (see also iBernardini et al.l 120091 ) . In order to check 
if this fit can be further refined, we started from the 
September 2007 observation, freezing the colder temper- 
ature at Tc — 0.15 keV, and letting the hot cap tem- 
perature Tfi free to vary. Since the cap area Ah is not 
known a priori, nor it can be treated as a free parameter 
in our minimization scheme, we tried several values of 
A/i, corresponding to one up to eight patches of our sur- 
face grid (this means that Ah is n/32 of the star surface, 
with n = 1, . . . , 8). The best-estimate emitting area was 
then taken as the one giving the lowest reduced x^ &!■ 
the fit in the different trials. We verified that in all cases 
the same value of the cap area produces the minimum 
X^ in all energy bands. Independent of the emitting area 
chosen, we always found for Th a value compatible with 
'- 0.15 keV for both the September and March 2007 ob- 
servations. 

One can then conclude that, for these two epochs, the 
entire star is radiating at the same temperature, or if 
a hot cap exists, its area is smaller than ~ 3 % of the 
star surface (the size of our surface grid resolution) . For 
these epochs we report in table [2] the values of the cold 
temperature obtained by using the single temperature 



scenario. We note that the temperature value is at the 
border of our grid of parameter values, so that, strictly 
speaking, it should be regarded as an upper limit on Tc- 
However, we verified that the x^ steeply grows when Tc 
increases above 0.15 — 0.16keV. Although there is no 
guarantee that the same is true when Tc decreases, in the 
following we assume that Tc ^ 0.15 keV is a satisfactory 
estimate for the uniform temperature at these epochs. 

We then proceeded backwards in time, from Septem- 
ber 2006 till September 2003. Again, the cooler tem- 
perature is kept fixed while several values of Ah are 
tried. However, to account for the possibility that 
also Tc varies, we repeated the calculation for Tc — 
0.15, 0.20, 0.25, 0.30 keV, looking for the pair {Tc, Ah) 
which gives the lowest x^- Results are summarized in 
table [21 Although the fits improve with respect to the 
one-temperature model (see table S]) , the two geometri- 
cal angles still change from one observation to another 
and also across different bands at the same epoch. 

In order to reproduce more accurately the star thermal 
map, we generated a new model archive, increasing the 
number of surface patches to Nq x Ng, = 50 x 4 = 200. 
The temperature, electron velocity and twist angle are 
in the range 0.15 keV < T < 0.9 keV (step 0.15 keV), 
0.1 < /3 < 0.9 (step 0.2) and 0.4 rad < A(J)n-s < 1.2 rad 
(step 0.2 rad), respectively. We then assumed that the 
star surface is divided into three zones: a hot cap at tem- 
perature Th, a concentric warm corona at T^, and the 
remaining part of the neutron star surface at a cooler 
temperature, Tc- Again, we began our analysis from the 
2007 observations, fixing Tc — 0.15 keV, and search- 
ing for the value of the warm temperature T^,. Every 
fit was repeated for twelve values of the emitting area 
A^ = 0.5%, 1%, 2%,4%, ..., 20% the total surface. 
We found that the reduced x^ improves with the addi- 
tion of a warm cap at T^ ~ 0.3 keV, accounting for 0.5 % 
of the neutron star surface (see table [S]). We stress that 
this value is below the resolution of our previous grid, so 
the two results are consistent with each other. 

We then considered the two 2006 observations; in the 
two-temperature model based on the previous archive, 
these were reasonably reproduced with Tc = 0.15 keV 
and Th ~ 0.3 keV (note that Th for the two-temperature 
model corresponds to T^ in the present case). For these 
two observations we repeated the fit, fixing Tc at 0.15 keV 
while leaving T^, free to vary. The size of the emit- 
ting area was estimated by following the same procedure 
discussed above. We found an almost constant value, 
T^ ~ 0.3 keV, between March 2006 and September 2007, 
while the emitting area decreases in time. Also, we found 
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TABLE 2 

XTE J1810-197 PARAMETERS AND THERMAL MAP (tWO-TEMPERATURE MODEL)^ 



Epoch 


Afliiv-s 


/3 


S(°) 


x(°) 


Th (keV) 


Ah{%) 


Tc (keV) 


Ae{%) 


Scp03 


0.70 ±0.01 


0.80 ±0.01 


22.7 ±0.5 


144.2 ± 0.6 


0.71 ±0.01 


25.0 ±3.1 


0.30 


75.0 ±3.1 


Sep04 


0.67 ±0.01 


0.62 ± 0.02 


20.7 ±1.0 


158.2 ±0.2 


0.55 ±0.01 


18.7 ±3.1 


0.30 


81.3 ±3.1 


Mar05 


0.61 ±0.01 


0.49 ±0.01 


21.6 ±0.4 


147.1 ±0.6 


0.67 ±0.01 


12.5 ±3.1 


0.25 


87.5 ±3.1 


Sep05 


0.47 ±0.01 


0.53 ±0.05 


23.0 ±0.1 


159.0 ±1.1 


0.42 ±0.01 


9.4 ±3.1 


0.25 


90.6 ±3.1 


Mar06 


0.49 ±0.01 


0.50 ±0.11 


23.5 ±0.4 


149.4 ±3.5 


0.28 ±0.01 


6.2 ±3.1 


0.15 


93.8 ±3.1 


Sep06 


0.43 ±0.01 


0.71 ±0.16 


21.4 ±0.3 


155.7 ±1.7 


0.28 ±0.01 


3.1 ±3.1 


0.15 


96.9 ±3.1 


Mar07 


0.45 ±0.01 


0.6 ±0.01 


29.8 ±0.1 


162.8 ±0.1 


_ 


_ 


0.16 ±0.01 


100.0 


Sep07 


0.48 ±0.01 


0.70 ± 0.08 


22.4 ±0.1 


163.0 ±1.7 


~ 


— 


0.15 ±0.01 


100.0 



^ Total energy band; parameters with no reported errors are fixed. Parameter errors are calculated by the min- 
imization routine mpcurvef it .pro, and are at la. Errors on the area correspond to the smallest patch of the 
grid. 



no need for a further component at Th at these epochs. 
On the other hand, results for the two-temperature case 
(see table [5]) show the presence of a component with 
temperature higher than 0.3 keV, in the period between 
September 2003 and September 2005 (while the cooler 
one varies between 0.25 and 0.30 keV). It is tempting to 
associate this to a transient hot cap that appears only in 
the first period after the outburst, superimposed to the 
other, longer-lived emitting zones. 

To test this possibility, we re-fitted the first four ob- 
servations by fixing the coldest temperature at Tc = 
0.15 keV, the warmer one at T^ = 0.3 keV, and leav- 
ing only the hotter temperature free to vary. For each 
observation the pulse profile fits were computed for ev- 
ery combination of Ah and Ayj chosen among the twelve 
values in the range 0.5 %- 20 % introduced before, and 
looking for the minimum of the reduced x^ . Results of 
the lightcurve fitting at different epochs are listed in ta- 
ble 13] and shown in Fig. SI while the reduced y^ for the 
three thermal distributions is reported in table |4l 

A worry may arise whether the best-fitting values ob- 
tained from the minimization routine correspond indeed 
to absolute minima of the reduced x^ . In order to check 
this, and visually inspect the shape of the y^ curve close 
to the solution, we computed and plotted the reduced x^ 
leaving, in turn, only one parameter free and freezing the 
remaining five at their best-fit values. This also allowed 
us to obtain a more reliable estimate of the parameter 
errors which were computed by looking, as usual, for the 
parameter change which corresponds to a Icr confidence 
level (and reported in table [3]). 

We found that all values obtained with the 
mpcurvef it .pro routine indeed correspond to minima 
of the reduced x^ curve, with the exception of the tem- 
perature(s), for which there are observations (or energy 
bands) with very flat x^ curves (see Fig. [5|). In par- 
ticular, for the September 2005 observation the curve 
obtained varying Th is flat in all the three energy bands. 
Also the curves relative to T^ for the September 2006, 
March 2007 and September 2007 observations have the 
same problem. This can be understood by noting that 
in all these observations the size of the hot /warm region 
accounts for only < 2 % of the total neutron star surface: 
temperature changes in such a small emitting area can 
hardly influence the fit. In addition, for the March 2005 
and March 2006 observations the reduced x^ curve rel- 
ative to one of the temperatures is flat, but this occurs 
only for one of the three energy bands. The first case 



concerns the hot temperature and the soft band, the sec- 
ond the warm temperature and the hard band. As we 
discussed above, when the hot (warm) area shrinks it 
affects little the pulse profile; this shows up first in the 
energy band in which its emission contributes less, i.e. 
the soft (hard) band. 

Given these findings, we concluded that lightcurve 
analysis by itself is unable to yield an unique temperature 
value for the September 2005, the September 2006, and 
both the 2007 observations. On the other hand, spectral 
analysis is more sensitive to temperature variations, so 
that it is possible to infer a temperature value also in 
these cases. As it will be discussed in the next section, 
by combining the two techniques we can remove most 
of the uncertainties and validate the three temperature 
model presented so far (see sec. 14.2.21 for details). 

There are several physical implications than can be 
drawn from our model. The TAXP is seen at an angle 
X = 148lg with respect to the spin axis. The mis- 
alignment between the spin axis and the magnetic axis 
is ^ = 2^°^^ . These values of the two angles, and the 
corresponding errors, are calculated from the weighted 
average in the three energy bands. To get a quantita- 
tive confirmation that x and ^ do not change in time, we 
fitted a constant through the values of each angle as de- 
rived from the lightcurves fitting at the different epochs 
and found that the null hypothesis probability is < 1 %. 
We note that, formally, the misalignment between the 
spin and the magnetic axis is compatible with being zero 
at the 3(7 level. Low values of ^ produce, however, mod- 
els with pulsed fractions quite smaller than the observed 
ones and, despite ^ ~ might be still statistically ac- 
ceptable, we regard this possibility as unlikely because 
the amplitude is the main feature which characterizes 
the pulse, as the PGA shows (see ij4.11 the first prin- 
cipal component, Zi, is, in fact, directly related to the 
amplitude) . 

It emerges a scenario in which, before the outburst, 
the NS surface radiates uniformly at a temperature Tc ^ 
0.15 keV. Soon after the burst the thermal map of XTE 
J1810-197 substantially changes. The region around the 
magnetic north pole is heated, reaches a temperature of 
~ 0.7 keV and covers an area ~ 8% of the total star sur- 
face. This hot spot is surrounded by a warmer corona at 
^ 0.3 keV, that covers a further ^ 16% of the surface. 
During the subsequent evolution, the hot cap decreases in 
size and temperature until the March 2006 observation, 
when it becomes too small and cold to be distinguished 
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TABLE 3 

XTE J1810-197 PARAMETERS AND THERMAL MAP (tHREE-TEMPERATURE MODEL)" 



Epoch 


A</)jv-s 


/9 


5(°) 


x(°) 


Th (keV) 


T» (keV) 


Tc (keV) 


Ah{%) 


A„(%) 


Sep 03 


0.80+nf 


o.7oj:°°» 


27.8+*-f 


145.3J:*;^ 


nK9+0.14 
U-O^_0.i4 


0.30 


0.15 


8. ±0.5 


16. ±0.5 


Sep 04 


0.79t°:g^ 


n 78+0.09 


16.21*:^ 


MO.sljg 


n 4Q+003 

U-4y_0.22 


0.30 


0.15 


6. ±0.5 


14. ±0.5 


Mar 05 


0.62t°:g^3 


0.5llH^ 


22.2t?3^ 


146.91^:9 


49+o°i 

'J-^^-0.04 


0.30 


0.15 


4. ±0.5 


14. ±0.5 


Sep 05 


r, c-q + 0.10 


0.50t°:\l 


21 4+l°-3 


154.4113^9 


0.52: 


0.30 


0.15 


2. ±0.5 


10. ±0.5 


Mar 06 


n ,f.+0.08 
0.46_(,o4 


0.73+°o^° 


18.5l«i° 2 


143.0+^-^ 


- 


o.29«:;^ 


0.15 


- 


6. ±0.5 


Sep 06 


n 54+004 


o.42lg:ii 


22 4+^2-^ 
^^•^-20.0 


150.2+1*53 


- 


0.301 


0.15 


- 


2. ±0.5 


Mar 07 


n 4Q+0-20 

'J-^^-o.os 


0.43l°:i^ 


of,Q+12.3 
■^'^•'^-20.0 


153.9ti°« 


- 


0.291 


0.15 


- 


0.5 ±0.5 


Sep 07 


„ 47+0.07 
U-4'_0.04 


0.50in^ 


99 7+16.4 
^^•'-20.0 


145.8116^* 


- 


0.311 


0.15 


- 


0.5 ±0.5 



^ Total energy band; parameters with no reported errors are fixed. Errors are computed from the x^ curve (see text 
for details) and are at la. No errors are reported when they could not be calculated (flat x^ curves) and errors on 
the area have the same meaning as in Tab. [2] 
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Fig. 4. — Parameters evolution for XTE J1810-197, three temperature model. Left (from top to bottom): twist angle (At/)), bulk velocity 
(/3) and area of the different emitting regions. Right: the two geometrical angles, x ^'^'^ C- Details as in fig. |3] 



TABLE 4 
Reduced x^ for XTE J1810-197 "^ 



Epoch 


^red 


v2 
^red 


^red 


v2 
^^red 


T 




(IT) 


(2T) 


(3T) 


(XSPEC) 


(keV) 


Sep 03 


1.72 


1.58 


0.12 


1.22 


- 


Sep 04 


0.66 


0.42 


0.36 


1.93 


- 


Mar 05 


1.02 


0.98 


0.79 


1.50 


- 


Sep 05 


1.06 


0.40 


0.39 


1.52 


n. c:q + 0.07 

U-0O_Q Qg 


Mar 06 


2.94 


1.70 


1.25 


1.34 




Sep 06 


0.94 


0.38 


0.35 


1.36 


n 01 +0.03 


Mar 07 


2.88 


2.88 


2.37 


1.08 


(, 20+0.04 


Sep 07 


1.12 


1.12 


0.96 


1.29 


•^1+001 



" First three columns: reduced x^ obtained from the 
lightcurves fitting for total the energy band (results for 
the other two bands are similar). Last two columns: 
reduced x^ obtained from the spectral fitting in XSPEC, 
and corresponding temperatures. The temperature was 
left free to vary only at those epochs and for those com- 
ponents for which the lightcurve analysis did not pro- 
duce an unique value. Errors for the temperature are 
at Icr. 



from the surrounding warm corona. The warm region 
remains almost constant until September 2005, then de- 
creases in size, and becomes a cap in March 2006, fol- 
lowing the hot spot disappearance. In September 2007 
(our last observation for XTE J1810-197) the warm cap 
is still visible, even if its area is down to only ~ 0.5% 
of the total. The twist angle is highest at the begin- 
ning of the outburst (September 2003) and then steadily 
decreases until it reaches a more or less constant value 
around September 2005. The electron velocity does not 
show large variations in time and stays about constant 
at /? ~ 0.5. 

Synthetic and observed lightcurves (in the total band) 
are shown in Fig [HI together with the fit residuals. We 
note that the residuals exhibit a well-defined, oscillatory 
pattern at all epochs. In our scenario, this can be pos- 
sibly associated to a more complicated thermal map, of 
which our 3T model is a first-order approximation (e.g. 
non-circular shape of the hotter regions, off-centering of 
the hot and warm areas). However, as discussed in some 
more detail in § [SJ no further refinement of the surface 
thermal map will be attempted here. Since XTE J1810- 
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Fig. 5. — Two examples of the different behavior of the reduced 
X^ for the warm temperature. The top curve refers to the March 
2006 observation (soft band) and exhibits a well-defined minimum. 
The bottom curve (September 2006 observation, total band) is so 
flat to make it impossible to gauge the best-fitting value and its 
errors. 
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Fig. 6. — Synthetic (3T model) and observed pulse profiles for 
XTE J1810-197 in the total energy band. Solid lines represent the 
best-fitting model, dots the observed lightcurves. Initial phases are 
arbitrary. The lower panel shows the residuals. 

197 pulse profiles are fairly sinusoidal, we can compute 
the pulsed fraction and its evolution in time at different 
energies. The comparison of model results with data is 
shown in Fig. [71 



Fig. 7.— The variation of XTE J1810-197 pulsed fraction with 
energy at different epochs. The red line refers to the model, blue 
dots to observations (errors are at Icr). 

4.2.2. Spectra 

In order to verify if the thermal map inferred from 
the pulse profile fits is reasonable, and in order to re- 
move the uncertainties in the value of the temperature 
at certain epochs (see § 14.2. ip . we examined the source 
spectra. The goal is to check if the parameters derived 
from our lightcurve analysis (twist angle, bulk velocity, 
size and temperature of the three emitting areas, and the 
angles Xj can also reproduce the spectral evolution 
of XTE J181 0-197 during the outburst decay. To this 
end, we use d the ntzanp; model that was implemented 
in XSPEC bv lNobih. Turolla. fc Zaiii ()2008l the model is 
not currently available in the public library, but it can be 
obtained from the authors upon request). The ntzang 
XSPEC model has the same free parameters as those used 
in our fits of the pulse profiles. In addition it contains the 
normalization and the column density. We caveat that, 
since this XSPEC model was created by assuming that the 
entire star surface emits at uniform temperature, strictly 
speaking is not directly suited to the present case. As an 
approximation, we fitted the spectra by adding together 
three (absorbed) ntzang models, each associated to one 
of the three thermal components, at temperatures T/j, 
Tu] and Tc, respectively. At each epoch the fit was per- 
formed by freezing A(/)jv_5, /?, T, x and ^ at the values 
derived from the fit of the lightcurve in the total energy 
band (see ij4.2.1[) . while the three model normalizations 
(which are related to the emitting areas) were left free 
to vary. We also required that the column density, Nh-, 
is the same for all the three spectral components and for 
all epochs. Since for the September 2005, the Septem- 
ber 2006 and both the 2007 observations the lightcurve 
analysis did not return an unique value for the hotter 
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temperature, we also left this parameter free to vary in 
these four observations. In ah these cases, we found that 
the fit converges to a value of the temperature close to 
the best-fitting value obtained from the lightcurve anal- 
ysis (see table HI). Moreover, the reduced x^ significantly 
worsens by varying the temperature, meaning that the 
spectra are much more sensitive to the presence of these 
components. Results are shown in Fig. |8l while the re- 
duced x^ for the fits at the various epochs are reported 
in Tab. IH The value of the column density is found to be 
Nh = (7.73 ± 0.50) X 10^1 cm~ ^ compatible at the l.Sa 
level with the one obtained by iBernardini et al.l (|2009f ) 
with the 3 BB model, Nh = (6.3 ± 0.5) x lO^i cm'^. 
We remark that, in assessing the goodness of the fits, 
only the normalizations of the three components (plus 
Nh) are free to vary; all the other model parameters are 
frozen at the best values obtained from the pulse profile 
analysis. Under these conditions, we regard the agree- 
ment of our model with observed spectra as quite satis- 
factory. We note that the presence of systematic residu- 
als at high energies (above 7-8 keV) may be hinted in the 
fits of the three ear lier obse r vation s (see Fig. [5]). As dis- 
cussed by Bcrnardini et al.l ()2009f ) they may be related 
to a harder spectral component which is however only 
marginally significant {3.2a confidence level) and quite 
unconstrained. Given that the high-energy residuals are 
comparable to (or smaller than) those of the 3 BB model 
used by Bcrnardini ct al. ( 2009) , we conclude that a hard 
tail is not significant also in our modelling and we did not 
attempt to include it in our fits. 

We checked how the reduced x^ for the spectral fit 
changes when the (frozen) parameters are varied within 
~ 2(7 from their best- fit value (as from the pulse fitting). 
This has been done changing one parameter at a time. 
We found that indeed the x^ increases quite smoothly 
in response to the change of each parameter, with the 
exception of x and £,. This is not surprising, since we 
knew already that the spectrum is not much sensitive to 
the geometry. We also tried a fit leaving all the param- 
eters free, apart from the two geometrical angles which 
were held fixed at their best-fit values. The fit returns 
parameter values which are the same, within the errors, 
as those derived from the pulse fitting and comparable 
values of the reduced x^, implying that the solution we 
presented is indeed a global x^ minimum. The same pro- 
cedure and the same conclusions hold also in the case of 
CXOU 164710.2-455216 (see S HX^ . 

4.3. CXOU J164710.2-455216 

Having verified that our model can provide a rea- 
sonable interpretation for the post-outburst timing and 
spectral evolution of TAXP prototype XTE J1810-197, 
we applied it to CXOU J164710.2-455216, the other tran- 
sient AXP for which a large enough number of XMM- 
Newton observations covering the outburst decay are 
available (see table[5]for details) . During September 2006 
the pn and MOS cameras were set in full window imaging 
mode with a thick filter (time resolution = 73.3 x 10~ 
s and 2.6 s for the pn and MOS, respectively), while 
all other observations were in a large and small window 
imaging mode with a medium filter (time resolution = 
4.76 X 10~^ s and 0.3 s for the pn and MOS, respec- 
tively). To extract more than 90% of the source counts. 





2 
Energy (keV) 






Fig. 8. — Spectral evolution in the eight XMM-Newton obser- 
vations of XTE J1810-197. Solid lines represent the model, while 
dotted lines refer to the single ntzang components (see text for 
details). Residuals are shown in the lower panel. 

we accumulated a one-dimensional image and fitted the 
ID photon distribution with a Gaussian. Then, we ex- 
tracted the source photons from a circular region of ra- 
dius 40" (smaller than the canonical 55", corresponding 
to 90% of the source photons, in order to minimize the 
contamination from nearby sources in the Westerlund 
1 cluster) centered at the Gaussian centroid. The back- 
ground for the spectral analysis was obtained (within the 
same pn CCD where the source lies and a different CCD 
for the MOS) from an annular region (inner and outer 
radii of 45"and 65", respectively) centered at the best 
source position. In the timing analysis, the background 
was estimated from a circular region of the same size as 
that of the source. EPIC-pn spectra were processed as 
in the case of XTE J1810-197 (see g21). 

4.3.1. Pulse profiles 

The analysis of the pulse profiles of C XOU J164710.2- 
455216 follows closely that presented in §4.2.11 In partic- 
ular, we first tried a single temperature and then a two- 
temperature model, encountering the same problems we 
found for XTE J1810-197. Finally, we applied a three- 
zone thermal map and this provided reasonable fits for 
the lightcurves, and the angles x and ^ were not found 
to vary in the same observation for the different energy 
bands and for different epochs. We did not attempt to 
fit the pulse profiles in the hard band after the February 
2007 observation because of the very few counts at ener- 
gies > 2 keV. As for XTE J1810-197 we started the anal- 
ysis from the last observation (August 2009) assuming a 
thermal map comprising a hot cap centered on the mag- 
netic pole at temperature Th, a concentric warm corona 
at T^ and the rest of the neutron star at the colder tem- 
perature Tc- Every fit was repeated for ten values of 
the hot cap area An = 0.5 %, 1 %, 2 %, 4 %, . . . , 16 % (of 
the total surface) and for 20 values of the warm corona 
area A^ = 0.5 %, 1 %, 2 %, 4 %, . . . , 30 %. Moreover, 
lightcurves fits were iterated for two values of the cold 
temperature Tc = 0.15,0.30 keV and also for two val- 
ues of the warm temperature T^, = 0.30, 0.45 keV. The 
hotter temperature was left free to vary. We found that 
in the last two observations, independent of the hot cap 
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TABLE 5 
CXOU J164710. 2-455216 XMM-Newton observations 



Label 


OBS ID 


Epoch 


Exposure time 


(s) 


total counts 


background counts 


Sep 06 


0311792001 


2006-09-22 


26780 




56934 


1709 


Feb 07 


0410580601 


2007-02-17 


14740 




18734 


1264 


Aug 07 


0505290201 


2007-08-19 


16020 




11710 


2384 


Feb 08 


0505290301 


2008-02-15 


9080 




4618 


1131 


Aug 08 


0555350101 


2008-08-20 


26360 




7357 


1689 


Aug 09 


0604380101 


2009-08-24 


33030 




4974 


1959 



"^ EPIC-pn 

size, Th is always ~ 0.45 keV, nearly indistinguishable 
from the temperature of the warm corona obtained from 
the fit. We concluded that, at least for our present sur- 
face grid resolution, in the last two observations there 
are only two thermal components that contribute to the 
emission, the cold and warm ones, and repeated the fit 
leaving T^, free to vary. Results are reported in table 
[S] and plotted in Fig. [HI while a comparison of the re- 
duced x^ for the three thermal distributions is given in 
table [71 Errors listed in the tables have the same mean- 
ing as in the case of XTE J1810-197. Again, when the 
spot at Th becomes very small its temperature can not 
be determined unambiguously. 

According to our model, CXOU J164710.2-4552116 is 
viewed at an angle x = 23^3 with respect to its spin 
axis. The spin and the magnetic axes are almost orthog- 
onal, £, = 84^i^^. This is a quite peculiar condition, and it 
seems to be the only one capable of explaining the char- 
acteristic three-peaked shape of the observed lightcurves 
within the present model. As for XTE J1810-197 values 
and errors for both angles are calculated as the weighted 
average of parameters in the three energy bands. Also in 
this case the probability that x and ^ are not constant 
in time is < 1%. 

Soon after the burst the thermal map of CXOU 
J164710. 2-455216 consists of three regions at different 
temperatures. The hottest region, around the north 
magnetic pole, has a temperature Th ~ 0.7 keV, and its 
area is ~ 8% of the total. This hot spot decreases in tem- 
perature and size as time elapses, until February 2008. 
In August 2008 the hot cap becomes so small in size and 
its temperature so close to that of the warm corona, that 
it is impossible to distinguish between the two regions. 
The warm corona has a temperature of ^ 0.45 keV, which 
remains about constant during the three years of obser- 
vations. In this case the corona area slightly increases 
with time, starting from ~ 20% and reaching ~ 30% of 
the NS surface. The third region has a lower temper- 
ature Tc ^ 0.15 keV and its area remains constant at 
~ 70% of the total. The twist angle is ^ 1.12 rad soon 
after the burst, and it decreases with time. There are 
hints that its decay is slower until August 2007, then 
proceeds faster. The electron velocity is about the same 
at all epoch {(3 ^ 0.2), apart from the last two observa- 
tions in which it strongly increases. This variation may 
be related to the change in the pulse shape (from three- 
peaked to single-peaked) and also to the increase of the 
pulsed fraction. 

A comparison between the observed and model pulse 
profiles is shown in Fig. 1101 Because of the inherent com- 
plexity and drastic time evolution of CXOU J164710.2- 



455216 lightcurves, the agreement is not as good as for 
XTE J1810-197. The fact that lightcurve fits return x^ 
values not much higher than those of XTE J181 0-197 
(compare Tab. |4] and [7|) refiects the larger uncertainties 
in the phase-binned source counts. We also note that the 
errors on the geometrical angles for CXOU J164710.2- 
455216 are smaller than those derived for XTE J1810- 
197, despite the worst agreement (see Tab. [3] and El). 
This is most probably due to the different shapes of the 
pulses in the two sources. Because of the very peculiar 
hghtcurve of CXOU J164710. 2-455216, which can be re- 
produced by our model only invoking a nearly orthogonal 
rotator, even small depatures of ^ and x from their best- 
fit values results in a rapid growth of the x^ . This does 
not occur for the rather sinusoidal pulse of XTE J1810- 
197 since the model can produce lightcurves of more or 
less the right shape in a wider range of angles. 

Besides being of limited use because of the complex 
shape of the pulse, the pulsed fraction analysis was hin- 
dered by the lower count rate, especially at low energies 
and was not pursued further for this source. As in XTE 
J1810-197, we checked that the values obtained from the 
minimization routine indeed correspond to minima of the 
reduced x^s. Again we froze five of the six parameters to 
the value obtained with the mpcurvef it .pro minimiza- 
tion routine, and calculated the reduced x^ around its 
minimum by varying the free parameter. The procedure 
was repeated for all parameters and all observations in 
the three energy bands. Again, for all parameters but the 
temperature, results obtained with the mpcurvef it .pro 
routine indeed correspond to the minima of the reduced 
X^ curve. There is one observation for which the x^ curve 
relative to the hot temperature is very fiat for all the en- 
ergy bands. This is the August 2008 observation, for 
which the size of the emitting area accounts for just 2% 
of the total neutron star surface. As in XTE J1810-197, 
we conclude that the fit is not very sensitive to the tem- 
perature variation for very small emitting areas. On the 
other hand, like in the previous case, it was possible to 
infer a value for the August 2008 hot temperature using 
the spectral analysis (see sect. I4.3.2|) . 

4.3.2. Spectra 

The spectral analysis for CXOU J164710.2-4552 16 was 
carried out using the same approach discussed in ^4.2.21 
We fitted three ntzang components, each representa- 
tive of an emitting region at different temperature, and 
froze all parameters apart from the three normalizations 
and Nh (which were forced to be the same for all the 
components and for all epochs). Moreover, since the 
lightcurve analysis the February 2008 observation failed 
to provide an unambiguous value for the hot tempera- 
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TABLE 6 
CXOU J164710. 2-455216 parameters and thermal map (three-temperature model)'' 



Epoch 


^<t>N~S 


P 


€(°) 


x(°) 


Th (keV) 


T„ (keV) 


Tc (keV) 


Ah(%) 


A„(%) 


Sep 06 


1 -IQ+O.OS 
^•^^-0.14 


n 1O-I-0.03 


83.5ti-« 


20.8«- 


0.70«:?? 


0.45 


0.15 


8. ±0.5 


22. ±0.5 


Feb 07 


1 07-f0.05 


1Q+002 


80.0+^-^ 


23.2+}-^ 


0.64tg:« 


0.45 


0.15 


6. ±0.5 


24. ±0.5 


Aug 07 


1 on+oi* 

^•'J'J-0.06 


lt;+002 


82.1+1-} 


19-2t?-3 


fiq+ois 

'J-O'^-0.10 


0.45 


0.15 


4. ±0.5 


26. ±0.5 


Feb 08 


n 77+0-21 

^•' '-0.12 


20+°''5 


85.7tll 


20.4t^:^7 


0.621 


0.45 


0.15 


2. ±0.5 


28. ±0.5 


Aug 08 


U.b5_oo7 


70+°"5 
'J- '^-0.05 


80.lt^-^ 


28.4+*-^ 


- 


49+''-°2 

u-^y-0.05 


0.15 


- 


30. ±0.5 


Aug 09 


0.55+_l:\l 


„ 7Q-I-0.06 


87.0t?o'9 


25.5tto' 


- 


46+''°5 


0.15 


- 


30. ±0.5 



" Total energy band. Errors have the same meaning as in Tab. [3] 
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Fig. 9. — Same as in Fig. [4]for the TAXP CXOU J164710. 2-455216; here time is computed starting from the September 2006 observation. 




TABLE 7 
Reduced x^ for CXOU ,1164710.2-455216^ 



Fig. 10.— Same as in Fig. [6]for CXOU J164710. 2-455216; com- 
puted pulse profiles refer to the 3T model and initial phases are 
arbitrary. 

ture, also this parameter was left free to vary. Results 
are shown in Fig. 1111 Given the approach we used for 
the fit, the agreement is quite satisfactory (reduced x^ 



Epoch 


v2 


^red 


^red 


v2 
Ared 


T 




(IT) 


(2T) 


(3T) 


(XSPEC) 


(keV) 


Sep 06 


1.05 


0.86 


0.31 


1.24 


- 


Feb 07 


1.32 


0.76 


0.65 


0.83 


- 


Aug 07 


0.97 


0.91 


0.44 


1.01 


- 


Feb 08 
Aug 08 


1.45 
1.45 


1.12 
1.23 


0.63 
0.79 


1.08 
1.23 


62+°°'^ 

U-0^_g Qg 


Aug 09 


2.03 


1.97 


1.52 


1.36 


- 



^ Same as in tab. |4] 

are listed if Tab. [7]). Systematic residuals at low (1- 
2 keV) energies are however present, especially in the 
September 2006, August 2007 and August 2008 obser- 
vations. Nh is found to be (2.14 ± O.OlSjj^iOi^ cm^, 
somewhat higher than that derived bv lNaik et al.l (|2008| ) . 
Nh = (1.73 ±0.03) x 10^^ cm^^. 

5. DISCUSSION 

The simultaneous study of the timing and spectral 
characteristics of the transient AXPs XTE J181 0-197 and 
CXOU J164710.2-455216 presented in this paper shows 
that the post-burst evolution of two sources share a num- 
ber of similar properties. In particular, the long-term 
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Energy (keV) 

Fig. 11. — Same as in Fig. |8]for the six XMM-Newton observa- 
tions of CXOU J164710.2-455216. 

variability of the pulse profiles and spectra appears to be 
(semi) quantitatively consistent with a scenario in which 
the star surface thermal distribution and magnetospheric 
properties progressively change in time. Our results 
were derived within the twisted magnetosphere model 
for magnetars and support a picture in which the twist 
affects only a small bundle of closed field lines around 
one of the magnetic poles. As discussed by Belob orodovl 
(2009), if the twist is initially confined along the mag- 
netic axis, the returning currents hit a limited portion of 
the star surface (typically a polar cap), which becomes 
hotter. In this scenario the post-outburst evolution is re- 
lated to the twist decay, during which the bundle shrinks, 
and the heated region decreases both in size and temper- 
ature. We found evidence for a cooling/shrinking of the 
heated polar cap in both sources, together with a de- 
crease of the twist angle. It should be noted that our 
magnetospheric model assumes a global twist, since no 
spectral calculations are currently available for a local- 
ized twist. 

Within this common framework, there are nonetheless 
differences between the two TAXPs. For XTE J1810-197 
we found that the star thermal map comprises three re- 
gions: a hot cap, a surrounding warm corona, and the 
rest of the surface at a colder temperature. The hot 
cap decreases in size and temperature until it becomes 
indistinguishable from the corona around March 2006. 
Also the warm corona shrinks, although its temperature 
stays about constant at ~ 0.3 keV. It becomes a cap in 
March 2006 and it is still visible in our last observation 
(September 2007) although its size is down to 0.5% of the 
entire surface. The rest of the surface remains at a tem- 
perature comparable to the quiescent one (as measured 
by ROSAT) during the entire evolution, indicating that 
the out burst likely involved on ly a fraction of the star 
surface. iBernardini et al.l (|2009[ ) obtained similar results 
using a 3BB model, although they did not attempt to 
locate the different emitting regions on the star surface 
nor to fit the pulse profiles. In their (purely spectral) 
analysis the hot region is visible slightly longer (until 
March 2006); the reason for the difference with respect 
to our results being most probably the resolution of our 
surface grid. Moreover, in our case the hot temperature 
decrease is more pronounced. The twist angle decreases 



from '^ 0.8 rad to '-^ 0.5 rad during the first two years, 
and then it remains roughly constant. 

Much as in the case of XTE J 18 10- 197, the thermal 
map of CXOU J164710. 2-455216 is well reproduced by 
three different regions. However, while the evolution of 
the hot cap is similar, i.e. it decreases in size and temper- 
ature until it disappears in the August 2008 observation, 
the behavior of the warm corona is different. Now the 
warm temperature remains constant at ^ 0.45 keV and 
the area increases. Actually, the area of the "hot-|-warm" 
region is constant and covers about ^ 30% of the surface, 
while the remaining '^ 70% is at a constant cooler tem- 
perature, ~ 0.15 keV. This is suggestive of a picture 
in which the "quiescent" state of the source is charac- 
terized by a two-temperature map, with a warm polar 
region superposed to the cooler surface. The outburst 
might have heated a portion of the warm cap, produc- 
ing the hot zone which then cooled off. It is intriguing 
to notice that the disappearance of the hot spot occurs 
at the same time (August 2008) at which the pulse pro- 
file dramatically changed, switching from a three-peaked 
to a single-peaked pattern. A quasi-sinusoidal shape of 
the light curve was observe d when the source was in qui- 
escence (jlsrael et al.l 120071 ). However, at that time the 
pulsed fraction was nearly 100% above 4 keV, likely in- 
dicating the presence of a small hot spot which is period- 
ically occulted as the star rotates. This is in agreement 
with our finding that this TAXP is a nearly orthogonal 
rotator. Whether CXOU J164710.2-455216 is presently 
approaching quiescence is unclear. If this is the case, its 
quiescent state is different from that observed in 2005 
and also from that of XTE J1810-197. 

It is worth stressing that our claim that the temper- 
ature does not change spatially in each of the regions 
should not be taken literally. The assumption that the 
surface can be divided in three (or two) thermal re- 
gions was mainly introduced to simplify the calculations 
while catching the essential features of the model. A 
smooth temperature variation within a zone is likely to 
be present. However, it is difficult to reconcile the ob- 
served pulsed fraction of XTE J1810-197 in the Septem- 
ber 2006 observation (> 10%, see Fig. ^ even accounting 
for the temperature gradient induced by the large-scale 
dipolar field. This may be an indication that, as our anal- 
ysis shows, there is a residual twist even in the quiescent 
state. 

In this respect we note that our spectral calculation 
is based on a rather fine subdivision of the star surface 
(50 X 4 patches in the final version of the archive), so 
we could have produced pulse profiles for arbitrary com- 
plicated thermal maps. The motivation of our choice of 
the thermal distribution (a hot polar cap and a warm a 
corona superimposed to the colder surface) is threefold: 
i) a model based on two thermal components, originating 
from a hot cap and a warm corona, was successfully ap- 
plied to XTE J1810-197 bv lPernafc Gotthelj (l2008l) : ii) 
inclusion of a third, colder component in the spectrum of 
the sa me source was shown to be statistically significant 
by B ernardini et al.l (|2009l ): and iii) it is consistent with 
theore tical predictions fo r a twisted magnetosphere in an 
AXP (|Beloborodovll2009l see above). In addition, this is 
the simplest map for which we were able to obtain con- 
stant values, to within the errors, for the two geometrical 
angles x ^^'^ ^ during the entire period covered by the 
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observations. 

In t heir analysis of XTE J1810-197. |Perna fc Gotthelj 
(|2008f) assumed that the X-rays come from two concen- 
tric regions with varying temperatures and areas, each 
emitting a blackbody spectrum; the rest of the surface 
was taken to be at zero temperature. They derived the 
angles x and ^, and, although their solution is not unique, 
they claim that the pair x ~ 53°, ^ ~ 23° is favored. 
While this value of ^ coincides with our estimate, the 
two values of the inclination of the line-of-sight are in 
substantial disagreement. Also the emitting areas of the 
hot/warm region and their temperatures turn out to be 
different in the two cases. Their estimate of the hot tem- 
perature is always higher than ours and the size of the 
warm corona is not monotonically decreasing. We re- 
mark that quantitative differences are to be expected 
given the different assumed sp ectral models (bla c kbody 
vs. RCS); moreover because iPerna fc GotthelB ()2008f ) 
did not include a colder region^. 

Finally, we caveat that our analysis relies on a num- 
ber of simplifying assumptions. We already mentioned 
that the synthetic spect ra we used were obta ined with 
the Monte Carlo code bv lNobili, Turolla, fc Z anc (20Q8), 
which was designed to solve radiation transport in a 
globally twisted magnetosphere. Even though we took 
thermal photons to originate mostly in a limited polar 
region, this does not self-consistently describe resonant 
up-scattering in a magnetosphere where only a limited 
bundle of field lin es is actually twist ed, as is probably 
the case in AXPs ( Beloborodovl [20091 ) . Moreover, as we 
discussed in ti4.2.21 the ntzEing XSPEC model is avail- 
able only in tabular form and it was created assuming 
emission at constant temperature from the entire star 
surface. As such, it is not suited to be applied directly 
to the present case. As a compromise, we decided to 
fit the spectra by adding together two/three (absorbed) 
ntzang components, each associated to one of the emit- 
ting regions, at temperatures T/j, T^ and Tc, respectively. 
While this procedure works (and is routinely employed) 
in the case of blackbody spectra, it is expected to be 
only approximately correct when different ntzang com- 
ponents are added together. The reason is that the ef- 
fects of resonant scattering on thermal photons depends 
on the location of the primary emission, since the mag- 
netospheric electron density is not isotropic. As a conse- 
quence, assuming thermal emission from a cap of limited 
size or from the entire star, even if the two are taken at 
the same temperature, will give rise to different spectra. 
We checked this approximation for all the spectra we an- 
alyzed, finding that the maximum relative error is ^ 0.6, 
while the energy-averaged error is always between 0.2 
and 0.4 both for XTE J1810-197 and CXOU J164710.2- 
455216. An example is shown in Fig[T3 Although we are 
aware that this is not optimal, it provides a reasonable 
way to describe radiation coming from a magnetar with 
non-uniform thermal emission within the context of our 
model. 




Fig. 12. — Comparison between the spectrum obtained adding 
three single ntzang model (red) and the spectrum of a neutron star 
with a thermal map consisting of three regions at different temper- 
atures (black). The two spectra are relative to the September 2004 
observation of XTE J1810-197. 

The monitoring of the two TAXPs XTE J1810-197 
and CXOU J164710. 2-455216, carried out with XMM- 
Newton in recent years, gave us the possibility to test 
the twisted magnetosphere model and understand how 
the physical parameters in the two sources change dur- 
ing the post-outburst evolution. We summarize our main 
findings below, remarking again that they were obtained 
under a number of assumptions (e.g. globally twisted 
field, three temperature thermal map). 

• Soon after the outburst onset the surface ther- 
mal distribution in XTE J1810-197 and CXOU 
J164710. 2-455216 is well described by three com- 
ponents: a hot cap, a surrounding warm corona 
while the rest of the neutron star surface is at a 
lower temperature. 

• The analysis of the pulse profile evolution for XTE 
J1810-197 revealed that both the hot cap and the 
warm corona decrease in size so that in the last 
observation (September 2007) virtually all the neu- 
tron star surface emits at a temperature compatible 
with the quiescent one. 

• The same analysis for CXOU J164710. 2-455216 
showed that the hot cap decreases in temperature 
and size, while the warm corona remains constant 
in temperature while it increases in size. In the last 
two observations we examined (August 2008 and 
August 2009) the source thermal map comprises a 
hot cap covering ^ 30% of the neutron star sur- 
face, while the remaining surface is cooler. There 
are hints that this could be the quiescent state of 
the TAXP. 

• For both sources the twist angle is highest at the 
outburst onset and then monotonically decreases 
in time until it reaches a nearly constant, non-zero 
value. 



6. CONCLUSIONS 



^ It was already noted bv IBernardini et al.l II2009I ) that the ad- 
dition of the colder component produces a monotonic decrease in 
both the hot and warm areas 



The same model configuration which best-fits the 
observed pulse profiles (thermal map, twist an- 
gle, electron bulk velocity, and geometrical angles) 
provides a reasonable description of XMM-Newton 
spectra in the 0.1-10 keV band for both sources. 
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To our knowledge this is the first time that a self- 
consistent spectral and timing analysis, based on a re- 
alistic modelling of resonant scattering, was carried out 
for magnetar sources, considering simultaneously a large 
number of datasets over a baseline of years. Present re- 
sults support to a picture in which only a limited portion 
of the magnetosphere was affected by the twist. Future 
developments will require detailed spectral calculations 
in a magnetosphere with a localized twist which decays 



in time. 



We are grateful to an anonymous referee for his/her 
constructive criticism and helpful suggestions which 
helped in improving a previous version of this paper. 
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AAE 1/088/06/0. 
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